Suppression of left-handed properties in disordered metamaterials 
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We study the effect of disorder on the effective magnetic response of composite left-handed meta- 
materials and their specific properties such as negative refraction. We show that relatively weak 
disorder in the split-ring resonators can reduce and even completely eliminate the frequency domain 
where the composite material demonstrates the left-handed properties. We introduce the concept 
of the order parameter to describe novel physics of this effect. 



Recently fabricated composite structures [J 0] °P en 
an unique possibility to study experimentally the prop- 
erties of the so-called left-handed metamaterials 
Specifically, the composite materials created by arrays 
of wires and split-ring resonators (SRRs) are known to 
possess a negative real part of the magnetic permeabil- 
ity and negative dielectric permittivity in the microwave 
range and they demonstrate many unusual prop- 

erties predicted theoretically long time ago 3, includ- 
ing negative refraction, inverse light pressure, reverse 
Doppler and Vavilov-Cherenkov effects, etc. Moreover, 
many ideas to create such left-handed metamaterials for 
optical frequencies are currently under study (see, e.g., 
Ref. @). 

The metamaterials created by arrays of metallic wires 
and microwave resonators are designed to operate for the 
wavelengths much larger then the period of the struc- 
ture, and their simultaneously negative dielectric per- 
mittivity and magnetic permeability occur in some fi- 
nite frequency range. Within the effective medium ap- 
proximation, dielectric permittivity becomes negative in 
a relatively large frequency domain due to a linear re- 
sponse of wires, whereas magnetic permeability of the 
structure can become negative in a relatively narrow 
frequency domain. It is this specific frequency domain 
where the induced magnetic momentum of each resonator 
is directed opposite to the external magnetic field being 
strong enough to produce collectively negative values of 
the effective magnetic permeability. 

In this Letter, we study the effect of weak disorder 
in the structure parameters on the magnetic response 
of composite and the existence of the frequency domain 
where the composite material exhibits the specific left- 
handed properties such as negative refraction. We show 
that even relatively small variation of the SRR param- 
eters can result in the dramatic decrease of the size of 
the left-handed frequency domain and, above a certain 
threshold value, disorder can suppress and even eliminate 
completely the left-handed properties of the metamate- 
rial. 

We study a composite structure made of metallic wires 
and single-ring microwave resonators, as shown schemat- 
ically in Fig.^a). This model is qualitatively similar to 
the model of double-ring resonators usually studied in the 
theory of left-handed composite media (see, e.g., Ref. 6] 



and reference therein). Within the effective medium ap- 
proximation, the dielectric permittivity of this structure 
can be calculated |7j, |8| and presented in the form 



eff(w) wl- 



LO (cj - ij e ) 



(1) 



where lo p — 27rc/dln (d/r w ), d is the size of the unit cell, 
r w is the radius of a metallic grid, j e = c /2aShi (d/r w ), 
a is conductivity of a wire, S is the effective wire cross- 
section, lo is the angular frequency, and c is the free-space 
speed of light. Usually, the effective plasma frequency, 
f p = lo p /2it, is between 10 GHz and 15 GHz. 

The most intriguing properties of the left-handed com- 
posite metamaterials are due to a SRR response which 
determines effective magnetic permeability. Magnetiza- 
tion of the metamaterial with three SRRs per unit cell, 
each of them having the axis in the perpendicular di- 
rection with respect to others as shown in Fig. [Ha), 
can be presented in the form (see, e.g., Refs. j^, lHj), 
M = x(cj)H', where 



r/LO 



V=l(^YlH8a/r) 



7/4] 



-1/2 



(2) 



(3) 



Here 
field, 



H' 



is 



(microscopic) magnetic 



[d g c 2 /ner 2 a(ln(8a/r) 



-7/4)] _ is 
a is a SRR radius, 
7/4] is damping coeffi- 



g is the size of the SRR slot, e is permittivity of 



acting 

LOQ = [u,g 

the SRR eigenfrequency, 
7m = c 2 /2a r S r [ln(8a/r) 
cient, d, 

the dielectric infilling the SRR slot in the structure 

Results (01 and © are obtained under the major as- 
sumption that all resonators in the structure are iden- 
tical. Now we consider the case when the size of the 
slot d g is random and it is characterized by a statistical 
distribution function. Our study can be motivated by a 
number of the recent efforts to create the metamaterials 
with the left-handed properties for shorter wavelengths 
and, for the smaller SRR size in the structure, the key 
parameters would be difficult to control in fabrication. 
In our structure, the parameter rj does not depend on 
d g [see Eq. ©], and only the SRR eigenfrequency lo is 
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affected by fluctuations of d g . Then, the equation for the 
magnetic susceptibility J2J) can be generalized to describe 
the case of randomly varying eigenfrequency, 



F(X) dX 



X 2 



(4) 



where F(X) is the normalized distribution function of 
the SRR eigenfrequencies, i.e. / °° F(X) dX = 1. In the 
standard case when all SRRs are identical, the distribu- 
tion function can be represented as F(X) = 5(X — uj ). 

To describe the coupling between the acting magnetic 
field H' and the macroscopic magnetic field H, we use 
the Lorentz-Lorenz formula [HI, H' = H + (47r/3)M, 
and present the effective magnetic permeability in terms 
of local magnetic susceptibility I0J as follows 



Ateff(w) 



1 + (8tt/3)xM 
1-(47t/3)xH" 



(5) 



For dcfinitcncss, we consider the Lorenz-type distribution 
of the SRR eigenfrequencies in the form 



F(X) 



(T/tt) 



(x - loo) 2 + r 2 ' 



(6) 



with a narrow width for r C wo such that the eigenfre- 
quencies of all resonators are close to some mean value uq. 
As a result, non-vanishing contributions to the integral 
are given by the values of X in the vicinity of X = u)q, 
and we introduce a new variable A, X = luq + A, where 
|A| <wo. We are interested in the behavior of the mag- 
netic susceptibility in the vicinity of to = loq, and we 
introduce to = loq + Q, where |fi| <C ujq. In this approxi- 
mation, Eq. |0J can be rewritten as 



xM 



2vr 



uj r 



dA 



^ (A 2 +r 2 )(A-ft + * 7m /2) 



(7) 



Using the contour integration, the integral in Eq. Q can 
be calculated explicitly, and the expression for the mag- 
netic susceptibility can be obtained in the following form: 



x(n) 



LOOT) 



1 



2 (-Q + iT)' 



(8) 



where f = V + j m /2. The result (JSJ shows that the ran- 
dom variation of the SRR eigenfrequencies is equivalent 
to some additional losses in the structure, and even for 
an ideal case of lossless resonators, i.e., when 7 m = 0, the 
composite structure with random variation of the SRR 
frequencies possesses effective losses. From the physical 
point of view, such effective losses resemble the collision- 
less Landau damping in plasmas |l2j| caused by the pres- 
ence of resonant particles and, simultaneously, it follows 
from the Kramers-Kronig relations. 

From Eqs. (J5J and JSJl, we find the frequency domain 
where the real part of magnetic permeability becomes 
negative, 0,% < D, < f2 2 , where 



ill- 



1. 



r c Tx/r 2 -r2, 



(9) 




FIG. 1: (a) Schematic of the metamaterial structure, (b) 
Dependence of the real part of the order parameter $lh on 
the effective combined losses F in the composite structure. 
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FIG. 2: Real part of the effective magnetic permeability as a 
function of the wave frequency for F/r c = 0.2 (solid), F/r c = 
0.5 (dashed), and f/r c = 1 (dotted). 



and r c = "kt\loq. Thus, the width of the frequency domain 
of the negative effective magnetic permeability is 



a^th = 2,/r 2 - f 2 



(10) 



The value T c has the meaning of the critical losses in 
the composite structure above which the domain with 
negative values of the magnetic permeability disappears, 
and the relative critical parameter T c /ujq depends only 
on the structure of the composite. Taking the character- 
istic values d = 0.5cm, a = 0.2cm, and r = 0.05cm, we 
obtain r c /cJo = 0.046. This result indicates that even 
for the case of lossless resonators, i.e. j m — 0, relative 
fluctuations of the SRR eigenfrequencies cannot exceed 
4.6%. Accordingly, if the slot size is d g = 0.01cm, (then 
loq/2tt = 9.6GHz) the critical variation can be estimated 
as Ad g ~ 5fim. Existence of such a critical distribution 
width places strict requirements on the manufacturing of 
such materials. Although for microwaves such require- 
ments can be easily met, they can result in substantial 
experimental difficulties for shorter wavelengths. In par- 
ticular, the recent proposal to fabricate left-handed meta- 
materials in optics using nanowires and 7r-shaped parti- 
cles 0, would require an accuracy better than 5 nm, 
that can be challenging for the existing technology. 

We note that premeditative introduction of disorder 
does not allow to increase the size of the left-handed fre- 
quency domain for any type of the distribution function. 
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Such a conclusion has a simple physical explanation. In- 
deed, disorder in the eigenfrequencies for any given SRR 
density results in a decrease of the effective number of 
resonators, which contribute into the negative magneti- 
zation. In the recent experiment |13| |. combining s-shapcd 
resonance particles, the authors fabricated the metama- 
terial with two domains of the negative refraction. Using 
the results obtained above, we may come to the con- 
clusion that the total size of the frequency domain with 
the negative magnetic permeability in the structure fab- 
ricated in Ref. [13| is less than it would be for the meta- 
material where all resonators are identical. Moreover, 
using the resonators with several different values of the 
eigenfrequency can eliminate completely the left-handed 
properties. 

Reduction and complete suppression of the frequency 
domain with negative magnetic permeability with a 
growth of the value of f can also be explained in a differ- 
ent way. Indeed, the left- and right-handed properties of 
the metamaterial can be treated as two different "phase 
states" of the structure. Thus, a transition from one state 
to the other can be interpreted as the phase transition 
of the second kind. The parameter T describes the effect 
of disorder in the structure, and it can be treated as "ef- 
fective temperature" . Such comparison is adequate since 
the effective temperature is determined as the width of 
statistical frequency fluctuations including both homoge- 
neous and inhomogeneous line broadening. To describe 
different phase states of the structure, we introduce the 
effective order parameter, as is usually done in the the- 
ory of phase transitions. In the absence of homogeneous 
(7m = 0) and inhomogeneous (T = 0) SRR line broaden- 
ing (i.e. in the absence of disorder), the metamaterial has 
the maximum width of the left-handed domain, which 
decreases and then disappears at the critical "effective 
temperature" T c termed as the temperature of the phase 
transition. We use the ratio of the left-handed frequency 
range to the maximum left-handed frequency domain as 
the order parameter, 

* LH = Afi LH /AQ^ ax) = Re[l - (f/r,) 2 ] 1 / 2 . (11) 

so that the metamaterial possesses the left-handed prop- 
erties when \1/lh ^ 0, i.e., below the effective critical tem- 
perature T c . Dependence of the order parameter "Jlh on 



f is shown in Fig. db) . The real part of the magnetic 
permeability as a function of the wave frequency is shown 
in Fig. [21 for different values of the ratio f/r c . For large 
values of r/r c , the frequency domain where the material 
possesses negative magnetic permeability is eliminated 
due to disorder. 

The result l|llf> suggests two ways for increasing 
the width of the effective frequency domain where the 
composite material possesses the left-handed properties, 
namely, (i) decreasing the effective temperature F, e.g. 
by improving the manufacturing technology for shorter 
wavelengths, and (ii) increasing the effective critical tem- 
perature T c , e.g. by a design of the resonators. 

Though the concept of the order parameter \&lh does 
not seem to be important in the context of the problem 
under study, it is expected to be useful in other (e.g. 
nonlinear) problems, when it becomes dependent on the 
field intensity, time, and coordinate. In this case, the 
order parameter will determine the material properties in 
different domains, describing the dynamics of the phase 
transitions between the states with left- and right-handed 
properties [g ■ The study of the nonlinear regime will be 
presented elsewhere. 

In conclusion, we have analyzed the effect of disor- 
der in the composite structures which exhibit left-handed 
properties. In particular, we have studied how random 
variation of the SRR eigenfrequencies in the structure 
can change the existence of the frequency domain where 
the effective dielectric permittivity and magnetic perme- 
ability are both negative. We have demonstrated that 
even relatively weak disorder in the SRR parameters can 
result in a dramatic reduction of the size of the left- 
handed frequency domain. More importantly, above a 
certain threshold value of this disorder, the left-handed 
frequency domain can disappear completely. We believe 
our result provide a useful guide for a design of novel 
types of metamaterials operating for shorter wavelengths 
where the effect of disorder is expected to be crucially 
important. 
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